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The long-wavelength properties of a noisy Kuramoto-Sivashinsky (KS) equation in 1+1 dimensions 
are investigated by use of the dynamic renormalization group (RG) and direct numerical simulations. 

Qlt is shown that the noisy KS equation is in the same universality class as the Kardar-Parisi-Zhang 
(KPZ) equation in the sense that they have scale invariant solutions with the same scaling exponents 
' in the long-wavelength limit. The RG analysis reveals that the RG flow for the parameters of the 

noisy KS equation rapidly approach the KPZ fixed point with increasing the strength of the noise. 
This is supplemented by the numerical simulations of the KS equation with a stochastic noise, in 
which the scaling behavior of the KPZ equation can be easily observed even in the moderate system 
size and time. 
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I. INTRODUCTION 



The most attractive models for surface roughening are the KPZ Q] and the KS 0,0 equations. The KPZ equation 
has a positive surface tension coefficient and is driven by a random forcing, while the KS equation is completely 
deterministic and is driven by inherent instabilities caused by a negative surface tension coefficient. In spite of such 
a difference between these equations, it was conjectured by Yakhot that the large-scale properties of the KS equation 
in 1+1 dimensions are described by the KPZ equation (a noisy Burgers equation) ^1,01, Ell- 
in order to clarify the claim, a number of numerical investigations in 1+1 dimensions have been developed 

The present understanding for this conjecture is that the spatiotemporal chaos generated by the negative surface 
tension becomes renormalized at long- wavelength into an effective positive surface tension and an effective noise term. 
In practice, the parameters in the effective stochastic equation that describes the long-wavelength properties of the 
KS equation in 1+1 dimensions were determined by the coarse-graining method 0,0 and other methods 0E|j an d 
on large scales it was shown that the KS equation behaves like the KPZ equation. However, Zaleski et al. were unable 
t-H ' to see the dynamic scaling phenomena which are numerically more costly to obtain precisely [3, |9|, [10( . Sneppen et al. 
observed the onset of crossover to asymptotic KPZ scaling by extensive numerical simulations on large system 0, 
but they could not find the KPZ dynamic scaling clearly. 

In the theoretical point of view, Yakhot used a RG approach [4| . Unfortunately, the theoretical assertion was based 
on uncertain ground, since the KS equation was treated with a perturbation theory around the unstable propagator 
which gives rise to an uncontrolled divergence. On the other hand, L'vov et al. proved the identity of the scaling 
behavior of the one-dimensional KPZ and KS equations in the long-wavelength limit under the locality of nonlinear 
interaction of these equations in the wave number space |ll| . They developed a perturbative treatment around the 
renormalized rather than the bare propagator by assuming the existence of such a renormalized propagator from the 
O ■ first. However, this approach needs to find a self-consistent scheme to examine its properties |l2j . 

There are other attempts to answer the claim addressed above, in which whether the KS equation with a stochastic 
noise term and the KPZ equation fall into the same universality class in 1+1 dimensions was investigated numerically 
and theoretically 0|. One is the numerical study of dynamic roughening in surfaces eroded by ion sputtering. 
The early and late time dynamics of an erosion model, which is inherently stochastic, were found to be the same as those 
obtained from the noisy KS equation |l3j . The other is the numerical simulation of a nonlinear stochastic equation 
describing the meandering of an isolated step on a crystal face grown from vapor. The nonlinear stochastic equation 
takes the form of the noisy KS equation above a critical supersaturation. The roughening function characterizing 
the step roughness obtained from the KS equation without a noise term coincides with that obtained from the one 
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with a noise term, with increasing distance from equilibrium (see FIG. 1 in [l4|)- This result suggests that the noisy 
KS equation exhibits the same behavior as the KS equation. Moreover, it is reported in Refs. 0, ^| that the 
steady-state spectrum (hkh-k) for the height h of a growing interface determined from the noisy KS equation was 
found to obey the generic k~ 2 scaling of the KPZ and deterministic KS equations for small wave number k, but the 
dynamic scaling properties were not clearly found in fl3| or not investigated in [T^ . This may be due to the smallness 
of the noise strength or limitation of the system size and time of the numerical simulations. In order to clarify the 
relation between the KS and KPZ equations from a different theoretical point of view, a dynamic RG analysis was 
performed for the KS equation with a nonconserved noise, and a stable fixed point of the RG flow equations for the 
parameters of the noisy KS equation was found, which was identified as the KPZ fixed point in Ref. |3- However, 
their values of the roughness and dynamic scaling exponents in 1+1 dimensions are different from those at the exact 
KPZ fixed point. 

The only difference between such a noisy version of the KS equation and the KS system is the effective noise term 
which originates from both deterministic noise (i.e., chaos) and stochastic noise [T^ |. It is important to understand 
how the interplay between these two noises determines the roughness of the surface. In this paper, we app ly the RG 
analysis to the KS equation with conserved and nonconserved noises to improve the RG results in Ref. [ija] . From the 
prediction of the RG results, we find a way to circumvent the limitation of the system size and time in the numerical 
simulations for the deterministic KS equation. 

The outline of this paper is as follows. In Sec. II, we perform the RG analysis for the noisy KS equation. In Sec. 
Ill, a part of the results predicted by the RG analysis is confirmed by the numerical simulations of the KS equation 
with a stochastic noise term. The conclusion is given in Sec. IV. 



II. RENORMALIZATION-GROUP ANALYSIS 



A. RG flow equations 



The noisy KS equation in one dimension is 

h t = vh xx - Kh xxxx + ^(h x ) 2 + r](x,t), (1) 

where the subscripts denote partial derivatives. Here h(x, t) describes the height profile of a one-dimensional surface 
above a substrate point x at time t, v is a negative surface tension coefficient, if is a positive surface diffusion 
coefficient, and A is the strength of the nonlinearity. 77(2;, t) is a Gaussian white noise with zero mean and the 
correlation 



(r](x,t)r]{x',t')) = 



d 2 

2D - 2D d — 

ox 



S(x-x')5(t-t'). (2) 



Here r\ is composed of the nonconserved and conserved noises whose strength is D and D^, respectively |l6j| . The 
Dd term in Eq. @ was not taken into account in Ref. ^jjj when performing the RG calculations. By introducing 
this term, we will obtain more reasonable RG results as stated in detail below. An equivalent equation to Eq. is 
obtained for u — —h x and / = — rj x as 

u t = vu xx - Ku xxxx - \uu x + /(x, t). (3) 

When A = 1, the nonlinear term in Eq. © is the same as that in the one-dimensional analog of the Navier-Stokes 
equation, then the variable u(x,t) can be interpreted as a one-dimensional velocity field in a compressible fluid 
The dynamic RG can be described through the Fourier modes with wave number k and frequency u>, in terms of 
which Eq. © takes the form 

{-iuj + isk 2 + Kk 4 )u{k 1 u J ) = f(k 1 u J )-—k[ ^ f —u{p,n)u{k- Pl uo-tt), (4) 

° ^bl<Ao -~ J -00 -~ 
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where Ao = tt/Ax is an upper cutoff in Fourier space, Ax = L/N is the lattice spacing in the real space, L and N 
being the system size and the number of grid, respectively. It follows from Eq. J2J that the Fourier transform of 
f(x,t) satisfies 

(f(k,oj)f(k\uj')) = 2{2n) 2 k 2 {D + D d k 2 )S(k + k')8{u + u'). (5) 
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First, following to the RG procedure 0, 0, El El 0] j we divide the velocity u(k, uj) into two components u > (fc, uj) 
and u < (fc,tt>), with the wave number satisfying A(Z) = A e~' < \k\ < A and |fc| < A(Z), I being a parameter. We 
eliminate (i.e., integrate away) the "fast" modes u > (fc,w), leading to an equation for the "slow" modes u < (fc,w) given 
by 



[-iuj + vk 2 +Kk 4 + Z(k,uj)}u<(k,uj) = f<(k,uj)-^k [ [°° ^ 

2 J\p\<MI) 2lT J-°° 2lT 
xu < (p,fl)u < (k — p,uj — CI), 

and f < (k,u>) satisfies 

(/< (k, uj)f< (k', (j)')) = 2(2vr) 2 fc 2 [D + D d k 2 + $(fc, u)] S(k + k')5(uj + uj'). 
E(fc, cj) in Eq. (j^J and <&(fc,w) in Eq. JJJ are given in the one-loop approximation as 

E(fc,w) = \ 2 k f % r ^[\G (p,Cl)\ 2 G (k-p,uj~Cl)p 2 (k~p)(D + D dP 2 ) 

■/A(0<|p|<A o Z7T J-oo 111 

+ \G (k - P ,oj- Cl)\ 2 G (p, n)(fc - pfp (D + D d {k ~ p) 2 )} , 



(6) 



(T) 



(8) 



and 



$(fc,w) 



dp 
2^ 



2^ 



Go(p,n)\ 2 \G (k-p,oj-Q)\ 2 



xp 2 (k-p) 2 (D + D dP 2 ) (D + D d (k-p) 2 ) 



(9) 



where Go(k,u>) — l/(—iu> + vk 2 + Kk A ) is the bare propagator. 

We are interested in the scaling behavior in the large system size and long time region, which correspond to k — > 
and uj — > limits, respectively. We first integrate the right hand sides of Eqs. (JBJ and over J7. Next, setting cj = 
and expanding the integrands in Eqs. JHJ and into Taylor's series up to order k A and k 2 , respectively, we can 
express them as E = 8vk 2 + 8Kk 4 and <f> = SD + SD d k 2 . Since the linear part of Eq. becomes unstable for v < 0, 
the bare propagator Go{k, uj) has a pole for the zero frequency at the wave number k — ko — (\v\jK) 1 ! 2 if A(Z) < feo- In 
order to avoid such an uncontrolled divergence, Procaccia et al. assumed the renormalized propagator and correlation 
function from the first, since there are rigorous proofs for the existence and boundedncss of the solutions of the KS 
equation 0,0]. Hence the propagator and correlation function in their E(fc,cj) and $(fc, u) are expressed in terms 
of the renormalized ones from the first. On the other hand in the RG method, if the integrations are performed over 
an infinitesimal wave number shell A (l — 51) < \k\ < A only, one can avoid the divergence due to the singularity 
of the bare propagator 15]. Then all calculations of 8v, SK, SD, and 5D d can be evaluated up to the first order 81 
without any problem, and we can define the renormalized parameters, v K = v + 8v , K < = K + 8K , D K = D + SD, 
= D d + 5D d . There is no correction to A, which is a consequence of the Galilean invariance y|, therefore, 

A. 

Second, in the RG procedure, we perform the rescaling, k = (1 + 8l)k, uj = (1 + z8l)u>, and u{k,uj) = [1 — 
(a + z)Sl]u < (k, lu), where a and z are the roughness and dynamic scaling exponents, respectively. The combined 
transformation yields the renormalized and rescaled parameters, v = [1 + (z — 2)8l\v < , K = [1 + {z — 4)Sl]K < , 
A = [1 + (a + z - 2)Sl}\<, D = [1 + (z - 2a - l)8l]D<, and D d = [1 + (z - 2a - 3)81}D<, where the variables with 
tilde denote the rescaled ones, while the variables without tilde are ones in the original scale. In the limit 81 — > 0, we 
obtain the one-loop RG flow equations describing the change in the parameters of the noisy KS equation under the 
RG transformation: 



and D< 

A< 



dv 
~dl 



G 



F(l + F) 3 



F 



(1-F) 



H 

G 



(10) 



dK 

~dT 



= K 



G 



2{l + Ff 



2Q — F + 2F l 



(2 — 2\F + 6F 2 + F 3 



H 

] G 



(11) 



dX 

H 



A [a 



2], 



(12) 
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dD 
~dT 



D 



z-2a-l + 



G 



{l + Ff 



H 
1+ G 



(13) 



dD d 
dl 



z - 2a - 3 



G 2 



2H(1 + F) 5 



x { 16 + 3F + F 2 + 2(9 - 5F)— + (2 - 13F - F 2 ) 



G 



(14) 



where we have defined the dimensionless coupling constants F(l) = v(l)/ (^K(l)Alj, G(l) = X(l) 2 D(l)/ (^Att K (I) 3 Aq^ , 

and H(l) = \{lfb d {l)/ [4nK(l) 3 Aiy which are expressed in terms of the rescaled variables. From Eqs. (|10|) - (|14fl . 
we can obtain the flow equations for F, G, and H: 

G 



dF 

dl 2(1 + F) 5 



TT 

6 - 12F + IIF 2 - F i + (2 + 19F 2 - 8F 3 - F 4 ) — 



(15) 



dG 
~dl 



— = 7G- 



G 2 



2(1 + F) 5 



TT 

76 - 7F + 4F 2 + 3F 3 + (2 - 71F + UF 2 + 3F 3 ) — 



-2(1 / (§ 



(16) 



G 2 



2(1 + F) 5 



1 16 + 3F + F 2 - (60 + 7F + 6F 2 + 3F 3 ) ^ 



(4-50F+19F 2 + 3F 3 



(17) 



There is another expression for the dimensionless coupling constants, /(/) = K (l)A 2 ] /v(l) 1 g(l) = 
X(l) 2 D(l)/ (47ri>(Z) 3 Ao), and h(l) = \(l) 2 Dd(l)A / (47r^(Z) 3 ) . However, they are not convenient when v is flowing 
towards zero, at which f(l), g(l), and h(l) diverge F(l), G(l), and H(l) can be expressed as F(l) = 1//(Z), 

G(l)=g(l)/f(l) 3 , andff(0 = M0//(0 3 - ' 



B. Fixed point of RG flow equations 



When putting £>d = (i.e., H = ) in Eqs. I|15|) and (|16fl . they do not reduce to the corresponding equations (20) 
and (19) at d = 1 in |15|. Therefore, we need to make some comments on the results in 15j, in which Cuerno et al. 
carried out the RG calculations for Eq. Q by taking into account only D term in Eq. (J2J, and found a stable fixed 
point at (F*,G*) = (-25.25,-722.8) at d = 1 with exponents z = 1.46, a = 0.54, and thus (3 = a/z = 0.37. They 
identified this as the KPZ fixed point. Their RG flow equations in the limit of \v\ <gC J^Aq reduce to Eq. (7) in [l7| . 
However, the negative values of F* and G* mean that unstable modes appear above the wave number k ~ {v* / \ K* \) 1 / 2 
because v* > 0, K* < at the fixed point. Moreover, the substitution of a — 0.54 in the scaling function of the 
energy spectrum yields E(k) cx /c~ - 08 for k — > 0, then we cannot say definitely that the energy spectrum for A; — > is 
independent of k (see the scaling solution of the e nerg y spectrum and the fluctuation-dissipation theorem below) . 

We have checked that the results at d = 1 in ,15| are obtained by making the substitutions p — * p + k/2 and 
fl — > Q + lj/2 in Eqs. 10 and ©, but without changing the integral region from A(l) < \p\ < Ao to A(l) < \p + k/2\ < 
Aq. We have performed the same integrations by replacing the integral region. As a result, we have obtained the 
same RG flow equations for v and D as Ref. but a different RG flow equation for K, since the effect of the 

replacement of the integral region appears when expanding the results up to order k 4 . We have carried out the same 
calculations of Eqs. JHJ and @ by directly integrating them without making the substitutions p — > p + k/2 and 
D, — > fi + u>/2, and confirmed to lead to the same results as those obtained by the replacement of both the variables 
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FIG. 1: The RG trajectories in the parameter space (F,G,H) for < I < 2.1, projected on (a) the (F,G) and (b) the (F,H) 
plane, for initial values v{Q) = -1, K(0) = 1, A = 1, D(0) = D (D = 0,0.1,5,40), and D d (0) = 0. 

and range of integration. The results are Eqs. 115JI and Ijlfijl without H term. Then we found a stable fixed point 
at (F*,G*) = (13.1868, 1064.43) with exponents z = 1.54, a = 0.46, and thus f3 = 0.30. We can recover the positive 
value of F* (i.e., K* > 0), therefore, it is found that the problem of the negative value of K at the fixed point is not 
attributed to the one-loop approximation |15| . Even if we made such a modification, the values of the exponents are 
still different from the exact values of the KPZ fixed point. 

In order to overcome this problem we have introduced Dd term in the noise correlator like Eq. (J2J). As a result, the 
parameter space (F, G) is extended to (F, G, H), and a stable fixed point for the RG flow equations (|T3|) - ((T7|l is found 
at {F*,G*,H*) = (10.7593,680.652,63.2614) with the scaling exponents, z = 1.5, a = 0.5, and thus = 1/3. Here 
the values of z and a are determined from equations dv jdl — and dD/dl — by the use of values of (F* , G* , H*). 
The scaling exponents are exactly those at the KPZ fixed point. Due to the Galilean invariance of Eq. A should 
remain unchanged under rescaling, therefore, Eq. I|12H leads to the scaling relation a + z = 2. Indeed, the values of z 
and a obtained above satisfy this scaling identity. 

Figure displays the trajectories of the dimensionless coupling constants (F(l),G(l), H(l)) under the RG trans- 
formation for an initial region < I < 2.1 with respect to four initial values (-1/Ag,0,0), (-1/Ag, 0.1/(4ttA£), 0) , 

(-1/A§,5/(4ttA£),0), and (-l/Ag,40/(47rAg),0), where we have set 9(0) = -1, K(0) = 1, A(0) = 1, D(0) = D 
(D = 0, 0.1, 5, 40), Dd(0) — 0, and Ao = ir/Ax — ir/0.5. These initial values will be used in the numerical simulations 
in the next section. Except for the case of (— 1/Aq,0, 0) (i.e., D(0) — Dd(0) = 0), the remaining three trajectories 
reach the same fixed point (F*,G*,H*) = (10.7593,680.652,63.2614). Moreover, it is found that the RG flow for 
(F(l),G(l), H(l)) rapidly approaches the KPZ fixed point with increasing the strength of D. The change in sign of 
F from F < to F > indicates that v changes from a negative value to a positive one. Once the value of D takes 
a different value from zero, the RG trajectory deviates from the F axis. For much smaller D the RG trajectories 
go along the F axis for much longer time (if we regard the parameter I as time), and then they turn the direction 
to the right and finally reach the KPZ fixed point. Therefore for very small D it takes much time (i.e., many RG 
transformations) to reach the KPZ fixed point. This behavior is also the case for D(0) = and F>d{0) ^ 0. We will 
numerically confirm in the next section that the KPZ scaling can be more easily observed by adding larger noises to 
the KS equation. In the case of D(Q) = Dd(Q) — 0, the initial value (— 1/Aq,0,0) goes to (— oo,0,0). If we see this 
behavior in the parameter space (f(l),g(l),h(l)), the initial value (— Aq,0, 0) approaches the fixed point (0,0,0). 

C. Undoing the rescaling 

Our real interest is in the effective theories without rescaling, so that, for comparison of the RG results with the nu- 
merical simulations in the next section, one must "undo" the rescalings (18| . In order to return to the original scale, we 
have to make the following transformations such that k = e~ l k, uj = e~ zl uj, v{l) = e _( ^ 2_2 ^i>(/), K(l) = e~( z ~ 4 '' K(l), 
X(l) = e- {a+z -^ r \{l), D(l) = e-^-^-^Dil), and D d {l) = e-^- 2a -^ l D d {l)- In terms of the variables in the original 
scale, the coupling constants can be expressed as F (I) = v(l)j (K(l)A(l) 2 ), G(l) = \(l) 2 D(l)/ {AirK(l) 3 K(l) 7 ), and 
H(l) = \(l) 2 D d {l)/ (4:TtK(1) 3 A(1) 5 ), where we have used A(Z) = A e~ l . Under the RG transformation the cutoff A is 
fixed, while the cutoff A(/) in the original scale corresponds to an appropriately chosen cutoff introduced in 0,0, in 
which short-wavelength degrees of freedom u(k) with |fc| > A(l) were explicitly eliminated. Going back to the original 
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FIG. 2: (a) The rescaled viscosity 0(1), (b) the rescaled noise strength D(l), (c) the viscosity v(l) in the original scale, and (d) 
the noise strength D(l) in the original scale, for 0(0) = v(0) = -1, K(0) = K(0) = 1, D(0) = D(0) = D = 0,0.1,5,40, and 
Dd(0) = D d (0) = 0. 



scale, the flow equations (TT3)l - <jT7|) remain unchanged, but the flow equations l|ll) |) -lfTl )l are changed to 



dv 



G ( ,H 

WTWV + f + (1 - f) g 



(18) 



dK 

~dJ 



= K 



1 26 - F + 2F 2 + F 3 + (2 — 21F + 6F 2 + F 3 )^ 



(19) 



rfA 
~dl 



0. 



(20) 



dD 

~dJ 



D 



G C HY 



(21) 



dD d 
dl 



Da 



G 2 



2H(1 + FY 



x J 16 + 3F + F 2 + 2(9 - 5F)^ + (2 - 13F - F 2 ) 



G 



G 



(22) 



These flow equations take the similar forms as Eqs. (I10 |I -(I14 [1 . but in the original scale it is not necessary to consider 
the contribution from the rescalings. 

Figures |21 (a) and (b) show the behavior of the rescaled viscosity 0(1) and noise strength D(l), which are obtained 
by solving Eqs. lfTU|l and lfH5|) by the use of Eqs. lfi"B|) - (fT?}l under the assumption that the values of a and z are given 
by a = 0.5 and z — 1.5. It is found that and D converge to steady-state values 0* and D* . Their values become 
large with increase of D(0) = D. On the other hand, Figs. |21 (c) and (d) show v(l) and D(l) in the original scale, 
which are obtained by solving Eqs. (|18fl and l|21|) by the use of Eqs. (|15fl - <|17|) or by making the scale transformation 
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v(l) = e~^ z ~ 2 ^ l 0(l) and D(l) — e~( z ~ 2a ~ 1 > l D(l). It is found that v(l) and D(l) increase with the parameter I, which 
indicates that the RG analysis cannot predict the values of the effective viscosity i/ c g and the effective noise strength 
D e g in the intermediate scaling region where the linear term dominates the nonlinearity in the effective equation 

obtained from Eq. © jlHIilia- B y 

using a discrete stochastic model of erosion processes by ion sputtering, it was 
found that a periodic ripple morphology characterizes the initial stage of the evolution, whereas the surface displays 
self-affine scaling in the later stage ^J. I n the RG procedure one must rescale the space, time, and height in order to 
obtain a system similar to the original, which is based on the assumption that the interface is self-affine. Therefore we 
can say that the RG analysis can determine the self-affine scaling exponents of rough interfaces, but cannot predict 
the crossover from the ripple structure to the rough surface. In Figs. |2] (a), (b), (c), and (d), we have chosen the 
initial values, 0(0) = v(0) = -1, K(0) = K(0) = 1, D(0) = D(0) = D = 0,0.1,5,40, D d (0) = D d (0) = 0, and 
A = 7r/ Ax = tt/0.5. It should be noted that the values of (F* , G* , H*) are universal in the sense that they do not 
depend on the initial values, while 0* and D* depend on them. 



D. Fluctuation-dissipation theorem 

From the linear part of the effective equation obtained from the noisy KS equation, the velocity correlation function 
in the original scale can be expressed as 0, 



(u(k,u;)u(k'^')) 2(D(l)k 2 + D d (l)k , 

(2-K)H(k + k>)5(u + u>) u 2 + (v(l)k 2 + K(l)k±f 
On the other hand, the velocity correlation function in terms of the rescaled variables is given by 



2 (b(l)k 2 +D d (l)k 4 ) 
C(k, Q) = — ^ ■ -L?. (24) 



+ [o(l)k 2 + K(l)k A 



C(k,ui) and C(k,u>) are related as C(k,uj) = e^ 1 ~ 2a ~ z ^ >l C(k,uj) by the scale transformation, then the scaling solution 
is given by C(k,u) — k x ~ 2a ~ z '^(ijJ /k z ), where ^(x) is a scaling function. Substituting the values z = 1.5 and a — 0.5, 
we obtain C(k,oj) = fc _1 ' 5 \& (oj/k 1 - 5 ). Thus, the scaling solution for the effective equation obtained from the noisy 
KS equation is the same as that for the KPZ equation [lj and the KS equation |IlLll9| in 1+1 dimensions. It follows 
from (|23|) that the energy spectrum of velocity in the original scale can be written as 



1 i ( k \ 2 

win f^rth \ D(l) + D d (l)k 2 D(l) ] + m_[ml (o^ 

E ( k ) = / 7T C fc ' w = n\ i is(i\i2 = 77 7K ; 7T> ( 25 ) 

J 2?r v(l) + K(l)k 2 v(l) 1,1/ k V 

1 + W) yw)) 

while the energy spectrum of velocity in terms of the rescaled variables is given by 

e(d = m±^m = M^fKiT. (26) 



0(l) + K(l)k 2 0(1) i , l ( % 

1 F{1) \A Q 



E(k) and E(k) are related as E(k) — e' 1 2a ^ l E(k) by the scale transformation, then the scaling solution is given 
by E(k) oc fc 1 " 2 " . At the hxed point, the values (F*,G*,H*) = (10.7593,680.652,63.2614) yield the estimates 
H*/G* w 0.093 and 1/F* w 0.093. Therefore E(k) at the fixed point is D*/0*. Using the values z = 1.5 and 
a = 0.5, we obtain v(l) = e 0-5(/-«o) £>* anc j = e o.5(l-l ) jj* ; w h er e we have used the fact that for l < I, 0(1) 

and D(l) take the steady-state values 0* and D* as shown in Figs. |U (a) and (b). Using the cutoff A(l) — Aoe - ', 
they can be expressed as v(l) = (A / A(l)) aM1 ' lo/l) 0* and D(l) = (A Q / A(l)) OM1 ~ lo/l) D*. Therefore in the limit 
I — > oo, i/(l) <~ A(/) -0 - 5 and D(l) ~ A(l)~°- 5 . This shows that the values of v and D in the original scale increase with 
decreasing the cutoff A(l) 0,123] . However, their ratio remains unchanged, D(l)/v(l) = D* jO* , which is found to be 
a scale independent quantity. Accordingly, E(k) also can be put as D* /0*, which indicates that in 1+1 dimensions 
the fluctuation-dissipation theorem holds in the noisy KS equation too in the long- wavelength limit. If we do not 
introduce D d term in the noise correlator, E(k) without H(l) in Eq. (|25|l depends on k because of 1/F(l) ^ 0. 
In the numerical simulations of the KS equation without noise term, the equal-time correlation function takes the 
steady-state value for k — > 0, S ■ The present RG analysis has shown that this is also the case in the noisy KS 
equation. 
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FIG. 3: Time evolution of Wit) for the KS equation without the noise term for L = 200, 000. Two straight dashed lines denote 
the line with exponent 0.25 and 0.29. 

III. NUMERICAL ANALYSIS 
A. Scaling exponents 

The long-wavelength properties of the KS equation are expected to behave like the KPZ equation. Sneppen et at 
performed a large-scale numerical simulation of the KS equation without the noise term Q. They investigated the 
dynamic scaling of ((h — (/i)) 2 ) 1 / 2 ~ t^ and found crossover from the Edward- Wilkinson scaling with (3 = 1/4. 
However, they could not find the KPZ scaling with (3 = 1/3 clearly. Much larger scale and longer time numerical 
simulations may be necessary to confirm the KPZ scaling. The RG analysis in the previous section shows that the 
fixed point (F* , G*, H*) is more easily attained by stochastic noises. We expect therefore that the KPZ scaling with 
(3 = 1/3 can be numerically observed in the noisy KS equation with moderate system size. 

In the usual numerical simulations of the KS equation without the noise term 0, la, m IliJ , all the parameters in Eg . 
© or © can be fixed to 1 by appropriate rescaling. The only control parameter is then the system size L |7ll8ll9l. Il0| . 
On the other hand, in the noisy version of the KS equation, D and Dd are added as the control parameters. We have 
performed direct numerical simulations of Eq. with v = — 1, K = 1, \ = 1, and investigated the time evolution 
of W{t) = ([h(x,t) - (/i(x,*))f) 1/2 ' We have investigated both cases of D ^ 0, D d = and D = 0, D d ^ 0, and 
obtained similar results. We will show only the results for D ^ and Dd = 0, since the case is more natural for the 
noisy surface growth 01 . We have used the Heun method for the numerical simulation with Ax = 0.5 and At — 0.005. 
(Sneppen et al. used the Euler method with Ax = 1 and At = 0.1, which might be too rough discretization.) The 
periodic boundary conditions are imposed, and the initial condition is h(x, 0) = 0. 

Figure displays the time evolution of W(t) for the KS equation for L = 200, 000. The ensemble average is taken 
for three runs. In this simulation, the noise term with D = 0.1 is added only for the initial time interval < t < 0.5, 
and the time evolution obeys the deterministic KS equation for t > 0.5. The double-logarithmic plot of W(t) shows 
that Wit) obeys the dynamic scaling t 1 / 4 fairly well. The slightly upward-curving line of W(t) probably represents 
the crossover to the KPZ scaling. The exponent near t ~ 50, 000 is about 0.29, which is still rather smaller than 1/3. 
That is, we could not confirm the KPZ scaling with (3 = 1/3 clearly, either. 

Figure 01 displays the time evolution of W(t) for the noisy KS equation with D = 0, 0.1, 5, 40 for L = 20,000. The 
ensemble average is taken for 10 runs. As the noise strength is increased, the time evolutions of W(t) are shifted 
upwards, since the fluctuations increase owing to the stochastic noises. For the deterministic case of D = 0, only 
the dynamic scaling with 1/4 is observed and even the crossover is hardly seen in this smaller size and shorter time 
simulation. For D = 0.1, the time evolution is almost the same as the case of D = 0. For D = 5, the exponent of the 
dynamic scaling increases from (3=1/4 towards (3 = 0.3. For D = 40, the exponent of the dynamic scaling is almost 
0.32, which is close to the exponent 1/3 of the KPZ scaling. This result is consistent with the RG result in Fig. ^ 
which shows that the KPZ fixed point is easily attained by the stochastic noises. If we can perform further extensive 
numerical simulation than that shown in Fig. the RG trajectory for the KS equation (i.e., the case of D = Dd = 0) 
would reach the KPZ fixed point. When D = Dd = 0, however, E in Eq. (JSJ and <3? in Eq. © are always zero. Then 
the RG trajectory is on the F axis in Fig ^ and never reach the KPZ fixed point. This is a qualitative difference 
between the RG analysis in the previous section and the numerical simulation. This is probably because the chaotic 
fluctuations generated by the deterministic KS equation are not well involved in the RG analysis when D = Dd = 0. 
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FIG. 4: Time evolutions of W{t) for the noisy KS equation at D = 0, 0.1, 5, and 40 for L = 20, 000. 

B. Modeling by the KPZ equation and estimate of the effective parameters 

Next, we analyze the numerical results of the noisy KS equation, based on the KPZ equation with the effective noise 
strength and the effective viscosity, assuming that the modeling by the effective KPZ equation is a good approximation 
for the behavior in the large spatial and temporal scales for the case of small D. The KPZ equation is written as 

h t = v c fth xx + ^(h x ) 2 + r](x,t), (27) 



(r](x, t)r}(x', t')) = 2D cS 5(x - x')6(t - «'). (28) 
Sneppen et al. evaluated the effective viscosity constant v c s and the effective noise strength D c g from the direct 

1 /2 

numerical simulations of the deterministic KS equation. They evaluated D c s/v c ^ using the relation 

2 = _p^i/ 2) (29) 

in the temporal range of the dynamic scaling with (3 = 1/4. In the long-wavelength region, the equilibriumlikc 
equipartition law for h x is observed for the KS equation 0, 0, S EH • Then, the second relation 

— =Lk 2 (\h(k,t)\ 2 ) (30) 

determines the ratio of D b r and u e g. From the two relations, Sneppen et al. evaluated -D c ff and f e g as D c g = 6.4 
and v e s = 10.5. They performed numerical simulation of the KS equation Q with v = — 1, K = 1, A = 2, and the 
above value of D e g is a rescaled value for A = 1. We note that the results by them are obtained by the replacements 
h — > 2h, W—> 2?7, and D c s — > 2D c g from the present results. The effective viscosity v e s was estimated with other 
methods @, El> and similar values of 6 < u e g < 10 were obtained. We have evaluated the effective viscosity v e s and 
the effective noise strength D e g for the noisy KS equation using the two relations l|29|l and i|30ll . We have performed 
numerical simulations of L = 20, 000 to evaluate the relation l|29|l . and numerical simulations of L — 4096 to evaluate 
the relation (|30|l . 

Figure[S](a) and (b) display v e g and D e s as a function of D. The effective viscosity v c $ decreases from v e $ ~ 8 near 
D = 0, and takes a value of about 5 for large D. The effective noise strength D e ff increases with D, since the stochastic 
noises are added to the deterministic noises generated spontaneously by the chaotic behavior. The values of v e g and 
D e s at D = are slightly different from the values estimated by Sneppen et al. Figure[S^c) displays the ratio D e ff/u e g 
as a function of D. The rhombi denote the numerical values determined by Eq. H30J) and the crosses denote the ratio 
D/v obtained by solving RG flow equations l|10|) and H13|) with Ao = tt/Ax — tt/0.5, assuming that the values of a 
and z are given by a = 0.5 and z — 1.5, and the initial values are £(0) = v = —1, K(Q) = K = 1, D(0) — D, and 
Dd{0) = Dd = 0. As we have shown in Fig|2(a) and (b), v and D approach steady-state values v* and D*. We have 
evaluated the ratio D/v at the steady-state. Even if returning to the original scale, the ratio remains unchanged, 
D e g/i> e g = D* jv* , as discussed in the previous section. It may be said that the RG analysis gives roughly approximate 
values of D c s/i> e ff for the long- wavelength fluctuations except for the range of small D. 
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FIG. 5: (a) Effective viscosity and (b) the effective noise strength numerically estimated using the two relations (1291 and 1301 1 as 
a function of D. (c) The ratio Deff/veii obtained by numerical simulations (denoted by rhombi) and the RG analysis (denoted 
by crosses). 



C. Approximation of the surface growth law using the RG analysis 

The increase of the effective noise strength is expected to facilitate the crossover to the KPZ scaling as we have 
shown in Fig. 0] We consider the dynamic scaling based on the KPZ equation. We assume that u e g and D e g are 
evaluated using Eqs. (|2T?)l and (|3*U|) . The RG flow equation for g(l) — \(l) 2 D cS (l)/(4:Tri> e f[(l) 3 A ) is given as 

d 9( l ) _ „m o„m2 



dl 



g(l)-2g(iy, (31) 



which can be obtained also from Eqs. (|10fl . i|12|) . and (|13fl by putting K = Dd = 0. Equation l|31|) is easily solved 
as g(l) = {1 + tanh(Z/2 — c)}/4, where c is determined from g(0) — A 2 D c g / (Airi/^Ao) = (1 — tanhc)/4 with the use 
of the initial values of j> c ff(0) = ^ c ft, D c g(0) = D c s, and A(0) = A. If we choose g(0) — 1/2, we obtain c = — oo. 
The substitution of c = — oo in the above exact solution of g(l) yields g(l) — 1/2 for any I. That is, if we choose the 
initial values of v e s, D c g, A, and Ao to satisfy g(0) = A 2 Z? c ff/(47rz^ ff Ao) = 1/2, the dynamic scaling with (3 = 1/3 
is easily observed even for small t, since the system is at the fixed point of the RG flow. The rescaled parameters 
i> e ff(l) and D e g(l) obey Eqs. (|lTj|) and (|T3|) by putting K = D d = 0. Returning to the original scale, the parameters 
v eS {l) = e- {z - 2)l v cS (l) and D eS (l) = e -(^^)i D eS (l) obey the equations 

dv eS {l) 

= vcft{l)g(l), 



dl 

dD oS {l) 
dl 



= D oS (l)g(l). (32) 



The effective viscosity v c g{l) and D e fi(l) increase as v e ${l) = v e s exp{/ g(l)dl} = v e s exp[{Z + 21ogcosh(Z/2 — c)}/4] 
and D c ff(l) = D e ff exp[{Z + 21ogcosh(//2 — c)}/4]. For I — + oo, the wave number k is scaled as k ~ Aoe _i , so 
I ~ log(Ao/fc). That is, the effective viscosity and noise strength in the original scale increase as v c s(k) ~ fc -1 / 2 and 
Dcs(k) ~ fc" 1 / 2 for small k. If we assume that v e s(l) and D e g(l) are the effective viscosity and the noise strength at 
the wave number k = Aoe - ', the time evolution of W(t) is expected to obey approximately 

W(t) 2 = Y{\h(k,t)\ 2 )~Y f e - 2 ^( fc ) fc2 ( t -*')2^ cff (fc)^ = Y Dc l {k \ ( 1 - e -*>Mk)k 2 t) (33) 
k k -* k VeS ( k > k V ' 

where the summation is taken for the discrete wave numbers k = k n = 2irn/ L (n = 1, • • • , N/2). We can evaluate the 
time evolution of W(t) using Eq. (|3"3)l . 

Figure[f)] (a) compares the numerical result of the KPZ equation l|27|l for L = 20,000, v e s = 10, and D e g = 5 with 
Eq. (|3"5jl for L — 20,000 and 1,000,000. Equation gives a fairly good estimate for W(t). The dynamic scaling 
with exponent (3 = 1/4 is seen even for the KPZ equation at ^ c ff = 10. The saturation effect of W(t) owing to the 
finiteness of L is also expressed by Eq. i|33|) . Figure (a) also shows that very large size and long time are necessary 
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FIG. 6: Comparisons of W(t) by direct numerical simulations (solid line) with Eq. ()3.'-ty (dashed lines) for (a) the KPZ equation 
with f e ff = 10, Deft = 5 and L = 20, 000, (b) the KS equation for L — 200, 000 and (c) the noisy KS equation for D — 5 and 
L = 20, 000. 

to see the crossover to the KPZ scaling for the large viscosity v e $ = 10. Figure 0(b) compares the numerical result 
of the KS equation at D = for L = 200, 000 (which is shown in Fig. [3| with Eq. (|33[) shown by the dashed line for 
L = 200,000 and v e g = 7.85, D c g = 7.65, where v c g and D e g are the values at D = in Fig. El (a) and (b). Very 
slow crossover towards the KPZ scaling can be predicted by Eq. . Figure (c) compares the numerical results of 
the noisy KS equation at D = 5 for L — 20, 000 (which is shown in Fig. with Eq. I|33|) shown by the dashed line for 
L = 20,000, v e g = 4.65, D c g = 26.7, where u e g and D e g are the values at D = 5 in Fig. 03(a) and (b). The crossover 
towards the KPZ equation is clearly seen even in the smaller-size system. 

The corresponding summation formula for W(t) using the RG flows of Eqs. I|15|) - (|22|l is written as 

W i t f ~ J2 f e- 2Mk)k2+K{k)ki}( - t ' t ' ) {2D(k) + 2D d (k)k 2 }dt' 
ft Jo 

= y, D(k) + D d (k)k 2 / 2Mk )k* + K<k)k*\t\ (M) 

Y v{k)k 2 +K(k)k i \ )' K ' 

The growth law with exponent (3 — 1/4 and the crossover towards the exponent (3 = 1/3 did not appear in this 
formula. We suspect that fluctuations by the deterministic chaos are well involved in the RG analysis of Eqs. (|1(J|) - 
(f 221) especially for small D. That is a reason why we have assumed the KPZ equation for the basis of the dynamic 
scaling for small D. It might be related to a fact that the effective viscosities around 5 ~ 10 did not appear in the 
RG flows. 

However, if D is sufficiently large, the stochastic noises dominate the chaotic fluctuations, and the RG flows of 
Eqs. (|15l) - (|22|l becomes more plausible. Besides, the effective viscosity and the effective noise strength cannot be 
evaluated for large D, since the intermediate region with the dynamic exponent (3 = 1/4 becomes invisible, and the 
analyses based on the KPZ equation cannot be applied. Here, we show the applicability of the RG flows for the noisy 
KS equation through Eq. I|34|l for large D. Figure compares the numerical results for L = 20, 000 and D — 40 shown 
in Fig. Q]with Eq. <|34|l . Fairly good agreement is seen for the large value of D = 40. 

IV. CONCLUSION 

We have performed the RG analysis for the KS equation with conserved and nonconserved noises in 1+1 dimensions. 
By introducing the conserved noise, first, we have obtained the values of the scaling exponents which are in good 
agreement with those at the KPZ fixed point. Second, we have shown that the fluctuation-dissipation theorem is 
exactly satisfied in the noisy KS equation in the long- wavelength limit. Therefore we can say that the long- wavelength 
properties of the noisy KS equation in 1+1 dimensions is described by the KPZ equation fairly well. However, the 
RG analysis for the noisy KS equation could not yet predict the effective viscosity and the noise strength in the 
intermediate scaling region. We have numerically evaluated the effective viscosity and the noise strength for the 
effective KPZ equation. We have shown that the KPZ scaling can be easily observed even in moderate-size numerical 
simulations of the KS equation under stochastic noises, owing to the increase of the effective noise strength. 

We have not succeeded in including the effect of the chaotic fluctuations into the RG flows and evaluating the 
effective viscosities for small D or D = 0. This problem also would arise in the other methods such as the direct 
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FIG. 7: Comparison of the time evolution of W(t) for the nosiy KS equation with D — 40 and L = 200,000 (solid line) with 
Eq. J3D (dashed line). 

interaction approximation (see Fig. 2a in 0) and the mode-coupling study of the KS equation, which were unable 
to obtain spectra with a hump. A rough estimate of the effective viscosity is v c s « I^u/t, wnere ^cejj is the size of the 
cellular structures in the KS system and t is the linear growth time of the most unstable mode [1,13 . It is recognized 
that the range of linealy unstable modes, which are seen in the prominent hump in the spectra, play a crucial role in 
the nonlinear dynamical behaivor of the KS equation |12|. Even if we are interested in the long- wavelength properties, 
we would not disregard this dynamically dominant range. In order to clarify whether the linearly unstable modes 
are incorporated by the RG method, it is worth to consider a model where the fourth-order derivative in Eqs. Q or 
© is replaced by a sixth-order one. In such a variant model, it is reported that the values of v c g and D e g become 
huge, but the fluctuation-dissipation theorem is still well satisfied 0,0. The huge value of v c g can be qualitatively 
estimated from u c s ~ ^oii/ 7 "- ^ ^ nc linearly unstable modes are well involved by the RG method, we would obtain the 
larger values of u e g and D c g for the variant model than those for the noisy KS equation. This is under investigation. 
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